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ABSTRACT 

Scaling relations between central black hole (BH) mass and host galaxy properties are of fundamental 
importance to studies of BH and galaxy evolution throughout cosmic time. Here we investigate the 
relationship between BH mass and host galaxy total stellar mass using a sample of 262 broad-line 
active galactic nuclei (AGN) in the nearby Universe {z < 0 . 055 ), as well as 79 galaxies with dynamical 
BH masses. The vast majority of our AGN sample is constructed using Sloan Digital Sky Survey 
spectroscopy and searching for Seyfert-like narrow-line ratios and broad Ha emission. BH masses are 
estimated using standard virial techniques. We also include a small number of dwarf galaxies with 
total stellar masses Msteiiar ^ 10®'^ Mq and a sub-sample of the reverberation-mapped AGNs. Total 
stellar masses of all 341 galaxies are calculated in the most consistent manner feasible using color- 
dependent mass-to-light ratios. We find a clear correlation between BH mass and total stellar mass for 
the AGN host galaxies, with Mbh oc Msteiian similar to that of early-type galaxies with dynamically- 
detected BHs. However, the relation defined by the AGNs has a normalization that is lower by more 
than an order of magnitude, with a BH-to-total stellar mass fraction of MBn/Afsteiiar 0 . 025 % across 
the stellar mass range 10 ® < Mgteiiar/Af© < 10 ^^. This result has significant implications for studies 
at high redshift and cosmological simulations in which stellar bulges cannot be resolved. 

Subject headings: galaxies: active - galaxies: evolution - galaxies: nuclei - galaxies: Seyfert 


1. INTRODUGTION 

A growing body of evidence suggests that su- 
permassive black hole (BH) masses scale with the 
large-scale properties of their host galaxies, pri¬ 
marily of the bulge component (e.g., bulge mass, 
velocit y dispersion, inf r ared l uminosity, see for in¬ 
stance iMaeorrian et al.l flQQSl: iGebhardt et all l2000al : 
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and references therein). These correlations on the one 
hand hint to a joint evolution of BHs and galaxies, and 
contain crucial information on the cosmic assembly of 
structures; on the other hand they provide a way to 
estimate BH masses via a proxy which is often much 
easier to measure. Extending these estimates to a 
statistical ensemble of galaxies allows one eventually to 
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In the local Universe, z = 0, benchmark BH 

masses are measured through direct methods, such 
as stellar and gas kinematics, and at the time of 
writing ^90 galaxies h a ve dynamical BH masse s (see 
iKormendv fc Hoi I20I31 IMcGonnell fc Mai 120131 and 
http: / /blackhole.berkeley.edu). The masses of the bulges 
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of the host galaxies can also be measured with good pre¬ 
cision, usi ng either photometric bulge/disk decomposi- 
tion (e.g.. IMarconi fc HuntI 1200311 or kinematical fitting 
(iHaring fc Rixl 120041) c oupled with assumptio ns of the 
mass-to-light ratio 1 see IKormendv fc Holl2013l for a dis¬ 
cussion) . 

Estimates of the relative mass of BHs and their host 
galaxies at higher redshift, which are of fundamental 
importance to establish the timing of their growth, do 
not have access to the same wealth of information avail¬ 
able locally. The BH masses are measured through in¬ 
direct methods, and their unc ertainties are discussed 
at length iii the li terature (e.g., IVestergaard fc PetersonI 
120061 : iShenI 1201311 . The host properties are also es¬ 
timated very d ifferently. Excep t for gravitationally 
lensed galaxies (IPeng et al.l 1200611 and/or HST images 
([Schramm fc SilvermanI 120131 : iPark et al.l l2015|l , which 
cannot provide a large statistical sample, reliable decom¬ 
position into the bulge and disk components are very 
difficult due to lack of spatia l resolution and sensitivity 
(but see iSanghvi et al1l2014D . Normally, the total stel¬ 
lar mass is used instead, estimated thro ugh assuming 
a mass-to-light ratio or SEP fitting le.g.. iJahnke et al 
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tryi ng to select quasars w ith bulge-dominated host galax¬ 
ies (iDecarli et al.l 1201011 to lessen the discrepancy be¬ 
tween bulge mass and total stellar mass. The scaling 
between BH and total stellar mass at high -2 is then of¬ 
ten compared to the scaling between BH and bulge mass 
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Figure 1. Spectrum of a broad-line AGN illustrating our fitting method. Top: the redshift-corrected spectrum is shown in black and the 
continuum plus absorption-line model is plotted in blue. Bottom: chunks of the emission-line spectrum (after subtracting the continuum 
and absorption-line model). The best-fit models for the emission line regions are shown in red. The individual narrow-line Gaussian 
components are plotted in yellow. Broad Hq and H/3 Gaussian components are plotted in dark blue. The residuals are shown in gray with 
a vertical offset for clarity. In the upper left-hand corner of the Ha and H/5 chunks we show the reduced values. The reduced x^ values 
from the fits not including a broad component are shown below in parenthesis for comparison. 


at z = 0. iLauer et alJ (|2007D expose important biases in¬ 
curred when investigating the potential evolution of BH 
scaling relations by comparing samples at different red- 
shifts with different selection criteria (e.g., AGN activity 
vs. host galaxy properties). 

In this paper we aim at quantifying on a local sam¬ 
ple the difference between using total stellar mass and 
bulge stellar mass to calculate the BH-to-host relation¬ 
ship. While the tightest correlation appears to be the 
one between the BH and the bulge, we wish to provide 
a benchmark for high-redshift studies which cannot avail 
themselves of bulge masses (or dynamical BH masses). 
We therefore investigate the relationship between BH 
mass and total stellar mass in a large sample of nearby 
(z ~ 0) broad-line AGNs using techniques for estimating 
BH and galaxy masses similar to those used at z > 0. 
Making use of active BHs has the added advantage of 
extending our sample to the lowest-mass BHs known i n 
galaxy nuclei (jReines et al.ll2013tlBaldassare et al.ll2015ll . 

2. SAMPLE OF BROAD-LINE AGN 

We construct our sample of broad-line AGNs by an¬ 
alyzing Sloan Digital Sky Survey (SDSS) spectra of 
~67,000 emission-line galaxies and searching for objects 
exhibiting broad Ha emission (signifying dense gas or¬ 
biting a massive BH) as well as narrow emission-line ra¬ 
tios indicating photoionization by an accreting massive 
BH. Our parent sample of emission-line galaxies is culled 
from the NASA-Sloan Atlas (NSA), which is based on 
the SDSS Data Re lease 8 (DR8) spectroscopic catalog 
(|Aihara et al.ll20ilT) . While we use the NS A for selecting 
our parent sample of galaxies, we use our own software 
to analyze the SDSS spectra and search for broad-line 
AGN. Distance estimates come from the zdist parame¬ 
ter in the NSA, which is based on t he SDSS NSA r edshif t 
and the peculiar velocity model of lWillick et al.l (1199711 . 
We assume Hq = 70 km s“^ Mpc“^. 

2.1. Parent Sample of Emission-line Galaxies 

The NSA catalog of nearby galaxies (z < 0.055) 
provides a reanalysis of SDSS optical photometry us¬ 


ing SDSS ugriz images with the i mproved background - 
subtraction technique described in iBlanton et al.r(l2011ll . 
The NSA also provides a reanalysis of spectroscopic 
data from the S PSS using the method s described in 
lYan fc BlantonI (12(31211 and lYanI (|20nll . We select 
emission-line galaxies in the NSA by imposing modest 
signal-to-noise (S/N) cuts on emission-line measurements 
reported in the NSA. We require S/N > 3 for the flux 
and S/N > 1 for the equivalent width (EW) of the Ha, 
[Nll]A6584, and [OlIl]A5007 emission lines. We also re¬ 
quire S/N > 2 for H/3 and the [Sll]AA6716,6731 doublet. 
This leaves us with a parent sample of 66,945 galaxies. 

2.2. Spectral Analysis and Selection of Broad-line 
AGNs 

We retrieved the SDSS spectra of our entire parent 
sample of galaxies and analyzed th em with customized 
software that is described in detail in IReines et al.l (|2013[1 
and briefly reviewed here for completeness. First we 
model and remove the stellar continuum and absorption 
lines from the host galaxy using simple stellar population 
model templates spanning a range of ages and metallic- 
ities. Next we model the narrow emission line profile 
based on the [Sll]AA6716,6731 doublet. Once we have a 
suitable model of the [Sll] doublet, we use it as a tem¬ 
plate for fitting the narrow emission lines in the Ha -I- 
[Nll]AA6548,6583 complex. We fit the Ha -I- [Nil] com¬ 
plex twice, first with the narrow lines only and then al¬ 
lowing a broad Ha component. We accept the fit with the 
broad Ha component if statistically warranted (reduced 

is improved by more than 50/^) and the FWHM of 
the broad Ha component is > 500 km s“^ after cor¬ 
recting for the fiber-dependent instrumental resolution 
(e.g., see Figure |T]). This FHWM requirement avoids 
severe contamination from intensely star forming galax¬ 
ies with moderately broadened bases on Ha. We also 

^ Reines et al. (2013) used a threshold of 20% since they were 
focused on dwarf galaxies with low-mass BHs, which can have weak 
broad Ha emission. Here we choose a higher threshold to help 
eliminate objects with marginally detected broad lines in higher 
mass galaxies. 
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Figure 2. Narrow- line diagnostic diagrams for sources with detectable broad Ha emission. We use the classification scheme outlined in 
IKewlev et al.l ()2QQ6lh Our sample of broad-line AGNs is restricted to objects with narrow line ratios placing them in both the AGN region 
of the [OIIIJ/H/3 vs. [NII]/Ha diagram and the Seyfert region of the [OIIIJ/H/3 vs. [SII]/Hq diagram (red points). The typical error for the 
red points is shown in the lower right corner of each plot. 


Table 1 

Sample of 244 Broad-line AGN 


ID 

(1) 

NSAID 

(2) 

SDSS Name 
(3) 

Plate-MJD-Fiber 

(4) 

2 :dist 

(5) 

Mi (Host) 
(6) 

g — i (Host) 

(7) 

log Mi, 
(8) 

log Mbh 

(9) 

1 

25955 

J000907.90-I-142755.8 

752-52251-320 

0.0422 

-21.81 

1.06 

10.68 

6.2 

2 

22075 

J004236.86-104922.0 

655-52162-58 

0.0424 

-21.10 

1.10 

10.43 

7.1 

3 

6452 

J012159.81-010224.3 

398-51789-10 

0.0548 

-22.77 

0.49 

10.47 

7.7 

4 

23318 

J021011.49-090335.5 

667-52163-506 

0.0419 

-22.60 

0.99 

10.92 

8.1 

5 

11183 

J024912.86-081525.7 

456-51910-77 

0.0296 

-20.64 

0.96 

10.11 

5.7 


Note. — Golumn 1: identification number assigned in this paper. Column 2: NSA identification number. Column 
3: SDSS name. Column 4: Plate-MJD-Fiber of analyzed spectra. Colu mn 5: zdist p arame ter in the NSA, which 
is based on the SDSS NSA redshift and the peculiar velocity model of IWillick et al.l 1)199711 . Column 6: absolute 
2 -band magnitude of the host galaxy. Colum n 7: Q — i color of t he host galaxy. Magnitudes and colors have been 
corrected for fore ground Galactic extinction ([Schlegel et al.l [199811 and the AGN contribution has been removed as 
described in ES Column 6: log host galaxy stellar mass in units of M 0 , corrected for AGN contribution. Uncertain¬ 
ties are on the order of 0.3 dex. Column 7: log black hole mass in units of Mq. Uncertainties are on the order of 0.5 dex. 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown 
here for guidance regarding its form and content.) 


measure fluxes of H/3, [Olll] A5007, and the [Sll] doublet 
to plac e objects on standard narrow-line diagnostic dia¬ 
grams (iBaldwin et al.lIlQSlHVeilleux fc Osterbrockl[T987l : 
IKewlev et al.ll200lL l200fiHKauffmann et al.ll2003[) . 

Figure [2] shows the [OlIl]/H/3 vs. [Nll]/Ha and 
[OlIl]/H/3 vs. [Sll]/Ha narrow-line diagnostic diagrams 
for all objects with detectable broad Ha emission. To 
minimize contamination from potential sources of broad 
Ha other than ionized gas orbiting a BH (e.g., super¬ 
novae in star-forming galaxies and shocks in LINERs), 
we restrict our sample of broad-line AGN to those sources 
falling in both the AGN region of the [OlIl]/H/3 vs. 
[Nil]/Ha diagram and the Seyfert region of the [OlIl]/H/3 
vs. [Sll]/Ha diagram. We also visually inspect each indi¬ 
vidual object and cut sources with poor spectral fits (e.g., 
due to complicated line profiles including double-peaked 
lines) that may lead to erroneous BH masses based on 
broad Ha (see below). A handful of objects are also ex¬ 
cluded for reasons described in tf2.4l leaving us with a 
final sample of 244 broad-line AGN (Tables |T] and [5]). 


2.3. Black Hole Masses and Luminosities 

Single-epoch spectroscopic BH m asses are routinely es¬ 
timated for broad-line AGN s fe.g.. iGreene fc He I2nn7bl : 
IVesterEfaard fc Osmed 120091 : iSchulze fc Wisotzld 1201011 . 
Under the assumption that the broad-line region (BLR) 
kinematics are dominated by the gravity of the BH, 
the BH mass is given by Mbh oc RAV'^/G. The 
average gas velocity is inferred from the width of a 
broad emission line (typically H/3) and the radius of the 
BLR is estimated from the radius-luminosity relation de¬ 
fined by reverbera t ion-m apped AGN (e.g., iKaspi et al.l 
120051 : iBentz et al.l 1201311 . The proportionality con¬ 
stant depends on the unknown geometry and orienta¬ 
tion of the BLR. While these par ameters have been 


seen to var" 

^ from object to object (IKoflatschnv 

1200,31: 

IBentz et ai.l 

20091:IDennev et al.ll2010l:lBarth et ai.l 

ml. 


a single scaling factor is generally adopted from cali¬ 
brating the ensemble of reverber ation-based BH masses 
to the Mrh ~ y* rel a tion fe.g.. iGebhardt et al.ll^OObl: 
iFerrarese et al.1 120011 : iNelson et al.l 12004 lOnken et al l 
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Table 2 

Sample of 244 Broad-line AGN: Emission Line Measurements 


ID 

{H/3)n 

(H/3)6 

[O III] 

[N II] 

(Ha)„ 

(Ha)6 

[N II] 

[S II] 

[S II] 

FWHM (Hq)6 




A5007 

A6548 



A6583 

A6716 

A6731 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

1 

245(13) 


1659(77) 

301(13) 

1192(31) 

757(63) 

892(39) 

250(16) 

219(11) 

1501 

2 

481(26) 

2343(147) 

4960(87) 

271(14) 

2022(103) 

11321(84) 

803(24) 

505(27) 

417(29) 

2121 

3 

664(47) 

2073(289) 

6841(203) 

892(60) 

3563(237) 

16872(168) 

2640(116) 

627(53) 

612(42) 

3281 

4 

494(33) 


1876(84) 

761(49) 

1802(116) 

6309(715) 

2252(150) 

748(57) 

585(69) 

7720 

5 

32(5) 

331(20) 

296(20) 

19(2) 

176(15) 

524(18) 

57(5) 

29(3) 

28(3) 

1081 


Note. — Column 1: identification number assigned in this paper. Columns 2-10: emission line fluxes with units of 10“^^ 
erg s“^ cm~^. Errors are given in parenthesis. We have not applied an extinction correction. The subscripts n and h indicate 
the narrow and broad components of the line, respectively. A three-dot ellipsis indicates that no line was detected. Column 11: 
FWHM of the broad Ha component. 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for 
guidance regarding its form and content.) 


2004 iGreene fc Hd 120061 : iPark et ^ l2012t iGrier et ^ 
mllHo fc Kimll2014n . 

We estimate BH masses for our sample of broad-line 
AGN using the single-epoch v i rial m ass estimator given 
by equation 5 in iReines et alJ (|2013D : 


log 




log e -I- 6.57 -I- 0.47 log 


-1-2.06 log 


^H„ 


1042 g,.g g-l 


/ FWHMhc 

\ 10^ km s“^ 


' 1 ) 


This e quation was derived following the approach out¬ 
lined in IGreene fc Hoi (j2005ll for using the broad Ha line, 
but incor porates the updated radius-luminosity relation¬ 
ship from iBentz et al.l (1201311 . Here we adopt e = 1.075, 
corresponding to t he mean virial factor < / >= 4.3 from 
Grier et ^ ( 201311 . where e = //4 (e.g., lOnken et al.l 
200411 . The distribution of BH masses for our sample 
of broad-line AGNs is shown in Figure [3 Viral BH mass 
estimates for broad-line AGNs are obvious ly very indi rect 
and carry uncertainties of ^0.5 dex (e.g.. lShenll2013ll . 

We estimate the bolometric luminosities of the AGNs 
using the conv ersion betwee n Lnrv and T 5100 given by 
Equation 1 in IGreene fc Hoi (120051 1. where I/Ha is the 
broad Ha luminosity and T5100 is the continuum lumi¬ 
nosity at 5100 A, and Lboi = IO. 3 L 5100 (|Richards et al.l 
I 2 OO 6 II . The range of bolometric luminosities is 41.5 < 
log Lboi ^ 44.4 and the median is log Lboi 43.4, ap¬ 
proximately 2.5 dex larger than the median of the distri¬ 
bution of broad Ha luminosities (see Figure [3]) • 


2.4. Total Stellar Masses of the Host Galaxies 

We estimate the total stellar masses, Mgteiiar, of galax¬ 
ies hosting broad-line AGN using mass-to-light ratios for 
i-band data (M/Lj) as a function ot g — i color following 
iZibetti et al.l (j2009l l. after removing the AGN contribu¬ 
tion to the integrated photometry. 

For each source, we estimate 5 -band and i-band flux 
densities of the AGN alone by constructing a mock AGN 
spectrum and convolving it with the SDSS filter through¬ 
put curves. The mock AGN spectrum consists of a power 
law (/a cx A“, where a = —1.56 for A < 5000 A and 
a = —0.45 for A > 5000 A; iVanden Berk et all 1200111 
plus the observed strong emission lines measured from 


the SDSS spectrum (H/?, [Olll], Ha, [Nil], [Sll]). We scale 
the mock AGN spectrum using th e conversion between 
L5100 and Lho (|Greene fc Holl200'^ . where the bro ad Ha 
luminosity is measured from the SDSS spectrum ( H2.2I) . 
For the vast majority of our sample (^85%), the AGN 
contribution to the total 5 -band and (-band flux densi¬ 
ties is less than 20% (see Figure |4]). To minimize erro¬ 
neous stellar mass estimates, we remove 5 sources from 
our sample in which the AGN dominates (i.e., contributes 
more than 50% to) the integrated flux densities. 

After removing the AGN contribution, the host-only 
AB magnitudes (corrected for Galactic reddening) are 
used to estimate galaxy stellar masses with a color - 
dependent mass-to-light ratio from IZibetti et ahl (|2009l l: 


log(M/L,) = 1.032(5 - () - 0.963 


( 2 ) 


We adopt a solar absolute (-band magnitude of 4.56 mag 
(|Bell et al.ll2003H . Errors on the stellar masses are ex¬ 
pected to be ^ 0.3 dex and are dominated b y uncertain¬ 
ties in stellar evolution (iConrov et al.ll200'9ll . 

In Figure|4|we compare total stellar masses of the host 
galaxies (with the AGN contribntion removed) to those 
derived from the integrated photometry (without remov¬ 
ing the AGN contribution). The median offset is 0.00 dex 
with a Ict scatter of 0.04 dex. In some cases, the stel¬ 
lar mass actually increases slightly since the 5 — ( color 
gets redder once the (blue) AGN is removed, and a red¬ 
der 5 — ( color increases M/Li. The effect of correcting 
for AGN contamination is minimal since our sample is 
dominated by Seyferts of modest luminosity ( )i2.3l) . The 
distribution of stellar masses is shown in Figure [S] 

We also compare the stellar m asses of our sample ba sed 
on the mass-to-light ratios in IZibetti e t al.1 (12 00911 t o 
those based on mass-to-light ratios in iBell et al.l ( ^03fl . 
as well as the stellar mass es provided in th e NS A based 
on SED fitting. For the IBell et aP (|2003ll masses, we 
again use M/Lt as a function of 5 — ( color. We then 
scale log M/Lj down by —0.093 dex (jZibetti et al.ll2009l : 
iGallazzi et al.l 120081 1 to account for the differences be- 
tween the Ghabrier initial mass function (IMF) used in 
and the scaled Salpeter IMF used 
The NS A stellar masses are de- 


iZibetti et al 


in IBell et ~ 


rived from the kcorrect code ([Blanton fc Row^l2007fl . 
which assumes a Ghabrier IMF and fits broadband op¬ 
tical fluxes from the SDSS and ultraviolet fluxes from 
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Figure 3. Distribution of the FWHM (left panel) and luminosity (middle panel) of broad Ha emission for our sample of nearby broad-line 
AGN. The distribution of virial BH masses calculated from equation 1 is shown in the right panel. 


the Galaxy Evolution Explorer when available. Stel¬ 
lar masses from the three methods (using integrated 
phot ometry in the N SA) are compared in Figure [51 
The [Bell et a h! (|2003[) m asses are systematically higher 
than the iZibetti et al.l (|2009l ) masses with a median 
offset of 0.21 dex {a = 0.10 dex) across the sample, 
with the largest discrepancies at low masses. The 
NSA (kcorr e ct) m asses are more consistent with the 
IZibetti et al.l (|2009l l masses with a median offset of 0.06 
dex (cr = 0.13 dex). 

3. ADDITIONAL OBJECTS 

We include the following additional objects in our in¬ 
vestigation of the relationship between BH mass and host 
galaxy total stellar mass. 

3.1. Dwarf Galaxies Hosting Broad-line AGN 

iReines et al.l (|2013l) carried out the first systematic 
search for AGNs in dwarf galaxies. The vast major¬ 
ity of objects in that sampl e are narrow- l ine A GNs 
and composites as defined bv iKewlev et al.l (|2006( ). for 
which we have no estimates of the BH masses. How¬ 
ever, 10 out of 136 AGN and composite objects have 
detectable broad Ha emission and narrow-line signa¬ 
tures suggesting the presence of an active BH (6 AGN 
and 4 composites). This subsample includes the well- 
studied d warf disk galaxy NGG 4395 hosting a Seyfert 
1 nuc leus ([Filippenko fc Sa,rgentlll989l: iFilippenko fc h 3 
120031) . two objects from the lGreene fc Hd l|2007al) sample 
of low-mass BHs. a nd the dwarf disk ga laxy presented in 
iDong et all (l2007| ). iB.eines et all (j2013D also provide in¬ 
formation on 15 galaxies with broad Ha emission, but 
narrow-line ratios dominated by star formation. As de¬ 
scribed in that work, there is likely significant contamina¬ 
tion from luminous Type H supernovae in that subsample 
and we therefore do not include those objects here. This 
issue will be further addressed in a forthcoming paper 
(V.F. Baldassare et ah, in preparation). BH masses and 
total stellar masses for the 10 br oad-line AGN and com¬ 
posites from IReines et al.l (|2013D are recomputec0 here 
in the same way as our full sample of broad-line AGN 
(see 112.31 and 112.41) and listed in Table |3l 

^ Here we usee = 1.075 in equation 1, rather than e = 1 as in 
IReines et al.l J201,'^ ). Stellar masses are computed us i ng eq uation 
[2] rather than taken from the NSA as in lR.eines et al.l 11201.11 ). 


We note that 3 of the broad-line objects from 
IReines et al.l (|2013D are also included in our main sample 
of broad-line AGN in this work (including NGG 4395). 
The remaining 7 objects are not recovered due to differ¬ 
ent selection criteria. In this work we exclude composite 
objects and also impose a more stringent threshold for 
flagging a source as having broad Ha emission (see 112.2p . 
since we are more concerned with having a clean sam¬ 
ple th an finding rare objects. In contrast, IReines et al.l 
(|20I3l) focused on finding low-mass BHs that can have 
weak broad Ha emission in dwarf galaxies that tend to 
have more active star formation. 

Follow-up high resolution spectroscopy of the 
IReines et al.l (|2013l) sample has led to the discovery 
of a new broad-line ob ject with Mbh 50,000 Mq 
([Baldassare et al.ll2015D . This object, designated RGG 
118 (object ID 118 in the Reines et al. paper), has the 
smallest BH reported in a galaxy nucleus. The BH 
mass for RGG 118 was estimated using equation 1 and 
here we estimate the stellar mass of RGG 118 to be 
ATsteiiar ~ 2.7 X 10® Mq using the SDSS i and g-band 
photometry in the NSA with equation |2l 

We also inc lude the well-studied dwarf Seyfert 1 
galaxy Pox 52 ([Barth et al.ll2004 iThornton et aill2008D . 
The mass of the BH in Pox 52 is Mbh 3 x 
10® Mq ([Thornton et al.l 120081) . As Pox 52 is not 
in the SDSS footprint, we estimate the galaxy stellar 
mass using H, Hand iA-band photo metry (correc t ed for 
the AGN contributio n ) pro vided bv I Bart hi et al.1 ([2004D 
and IThornton et al.l ([20081 ) with the following color - 
dependent mass-to-light ratio from IZibetti et ahl ([20091 ) : 

log(M/Lif) = 1.176(H-H) - 1.390. (3) 

We adopt a solar absolute AT-band magnitude of 3.32 mag 
([Bell et al.l 120031) . The resulting stellar mass for Pox 52 
is Mgteiiar 4.3 X 10® Mq. This is ^^2.8 times smaller 
than the stellar mass estimate d bvlTho rnton et al.l (120081) 
using the kco rrect code o f [Blant on fc BoweisI ([20071 ). 
We adopt the IZibetti et al.l (|2009[ ) mass for consistency 
with the rest of our sample. 

As with our main sample of broad-line AGNs ([[1|), we 
adopt uncertainties of 0.3 dex in stellar mass and 0.5 dex 
in BH mass for all of the dwarf galaxies discussed above. 
However we caution that the virial BH masses derived for 
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Figure 4. Top: Distributions of the ratio of AGN flux density 
to total (Host+AGN) flux density in the SDSS g and i bands for 
our sample of broad-line AGNs. Middle: Distributions of absolute 
p'-band magnitudes for our sample of broad-line AGN. The blue his¬ 
togram shows AGN-only magnitudes, the orange histogram shows 
host galaxy-only magnitudes and the black histogram shows the 
total (Host-hAGN) magnitudes. Bottom: Total stellar mass with¬ 
out correcting for AGN contamination vers us total stellar mass 
corrected for AGN contamination. See 3231 for details. The line 
shows the one-to-one relation. 



7 8 9 10 11 12 

log 


Figure 5. The distribution of host galaxy total stellar masses for 
our sample of broad-line AGN (corrected for AGN contamination) 
is shown in red. Our parent sample of emission-line galaxies is 
shown in blue, normalized to the number of galaxies in the red 
histogram (244 objects). We also show the mass distribution for the 
full NS A catalog (no emission-line cuts), again normalized by the 
number of galaxies in the red histogram. All masses were derived 
using g and i-band data in the NSA with the color-dependent mass- 
to-light ratio given in equation [2] 

these objects are based on an extrapolation from more 
massive and luminous AGNs and may carry additional 
err ors. We do not inc l ude the dwarf galaxies Henize 2- 
10 (iR.eine s et al.l 1 20111 : iR.eines fc Delleill2012t). Mrk 709 
(jR.eines et al.ll2014i) or ,11329-1-3234 (iSecrest et al.ll2015[) 
since the BH masses in these systems are uncertain by 
at least an order of magnitude. 

It is worth noting here that optical searches for AGN in 
dwarf galaxies suffer from severe selection effects. First, 
low-metallicity AGN, which are likely to reside in low- 
mass galaxies, have line ratios significantly different from 
Seyferts in more metal-rich systems and overlap with 
low-metallicity starbu rsts in the [OlIl]/H / 3 vs. [Nil]/Ha 
(i.e., BPT) diagrain (iGroves et al.l 120061 : iKewlev et alJ 
120131 : iReines et alJ 1201 41) . Moreover, broad Ha from 
small accreting BHs can be very weak and difficult to 
detect. Given the sensitivity of the SDSS s pectra and 
our search volume (z < 0.0551. iB.eines et HI (1201, 3[ ) esti¬ 
mate a minimum detectable BH mass of Mbh 10^ Mq 
if the BH is radiating at ^10% of its Eddington limit. 

3.2. Reverberation-Mapped AGN 

The most reliable A GN BH masses come from rever- 
berat i on mapping (e . g..lPeters on et al.l l2004 iBentz et al.l 
120091: iDenney et al.l l^lOl : iBarth et al.l l2011[) . Deter- 
mining the time lag between the continuum flux and 
broad emission line variability gives the light travel 
time across the BLR, and in turn the BLR radius 
when multiplied by the speed of light. The BLR 
radius-luminosity correlation derived from the sample of 
reverberation-mapped AGN (^50 objects) makes single¬ 
epoch virial BH masses for AGNs, such as th e ones 
used in this wor k, possible (e.g., iKasoi et~aI1 120051 : 
IBentz et al.l 1201311 . The reverberation-mapped AGN 
also provide a link between single-epoch spectroscopic 
BH masses and dynamical BH masses, as the ensem- 
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Table 3 

Additional Objects 


Name 

log M^, 

log Mbh 

(1) 

(2) 

(3) 


Dwarf Galaxies with Broad-line AGN 

RGG 1 

9.30 

5.44(0.50) 

RGG 9 

9.24 

5.00(0.50) 

RGG 11 

9.30 

5.29(0.50) 

RGG 20 

9.29 

6.10(0.50) 

RGG 21 

9.45 

5.80(0.50) 

RGG 32 

8.90 

5.28(0.50) 

RGG 48 

8.96 

5.18(0.50) 

RGG 118^^ 

9.43 

4.70(0.50) 

RGG 119 

9.12 

5.42(0.50) 

RGG 123 

9.36 

5.79(0.50) 

RGG 127 

9.36 

5.21(0.50) 

Pox 52^ 

8.63 

5.48(0.50) 

Reverberation-Mapped AGN 

Mrk 590 

11.38 

7.57(0.07) 

Mrk 79 

10.41 

7.61(0.12) 

Mrk 110 

10.14 

7.29(0.10) 

NGG 3227 

10.36 

6.77(0.10) 

SBS 1116-I-583A 

10.00 

6.56(0.08) 

Arp 151 

9.90 

6.67(0.05) 

Mrk 1310 

9.62 

6.21(0.08) 

NGG 4051 

10.16 

6.13(0.14) 

Mrk 202 

9.92 

6.13(0.17) 

NGG 4253 

10.34 

6.82(0.05) 

NGG 4395 

8.90 

5.45(0.14) 

NGG 5273 

10.25 

6.66(0.16) 

NGG 5548 

10.79 

7.72(0.02) 

Mrk 817 

9.87 

7.59(0.07) 

Mrk 290 

9.52 

7.28(0.06) 


Galaxies with Dynamical BH Masses 


M 32 

8.77 

6.39(0.18) 

NGG 1316 

11.48 

8.23(0.07) 

NGG 1332 

10.92 

9.17(0.06) 

NGG 1374 

10.33 

8.77(0.04) 

NGG 1399 

11.17 

8.94(0.33) 

NGG 1407 

11.43 

9.67(0.05) 

NGG 1550 

11.02 

9.59(0.07) 

NGG 2960 

10.72 

7.03(0.02) 

NGG 3091 

11.29 

9.57(0.04) 

NGG 3377 

10.14 

8.25(0.23) 


Note. — Column 1: object name. 
Dwarf galaxies with the des i gnatio n 
RGG are from IReines et al.l l|2013l j. 
Column 2: log host galaxy stellar mass 
in units of Mq. The AGN contribution 
has been removed for the dwarf galaxies 
and reverberation-mapped AGN. All 
stellar masses are estimated using color- 
dependent mass-to -light ratios from 
IZibetti et al.l (|20Q9l h Uncertainties are 
on the order of 0.3 dex. Column 3: log 
black hole mass in units of Mq. BH 
masses for the reverberation-mapped 
AGN are taken from th e AGN Black 
Hole Mass Database JBentz & Katzl 
[201511 . Dynamical BH masse s are taken 
from IKormendv fc Hd ll2Q13lh 

(This table is available in its en¬ 
tirety in a machine-readable form in 
the online journal. A portion is shown 
here for guidance regarding its form and 
content.) 

^ BH mass from Baldassare et al. (2015). 
^ BH mass from IThornton et al.l 11200811 . 



log (M 3 teiia/M 5 ), Zibetti+09 



log (M 3 teiia/MQ), Zibetti+09 


Figure 6. Comparison of stellar masses of our sample of broad¬ 
line AGN derived from different methods. We calculate mass-to- 
light ratios fo r the SPSS 2 -band d ata as a function of q — i color 
following both IZibetti et aLl (1200911 and IBell et al.l ll2003l'l and using 
integrated photometry provided in the NSA (AGN contribution 
has not been removed). NSA masse s are derived from the kco rrect 
code, which is described in detail in IBlanton &: Rowei^ d200'/T) . The 
lines show the one-to-one relation. 


ble of reverberation-mapped B H masses is calibrated 
to the Mrh — y* rel a tion fe.g.. iGebbardt et al.ll20d0bl: 


Ferrarese et al.l 120011: [Nelson et al.[ [2004 [Onken et al. 

to 

0 

0 

: Greene & Hoi 20061 [Park et al.l 120121 [Grier et al. 

2013 

ilHo & Kim[l2m 4[). 


In this work we include 15 reverberation-mapped 
AGNs with BH masses provided in th e AGN Black Hole 
Mass Database (jBentz &: Katd 1201511 . We adopt BH 
mas ses calculated wi th a mean virial factor of < / >= 
4.3 (IGrier et al.ll2013il as we did for the single-epoch spec¬ 
troscopic masses given by equation 1. We first cross¬ 
match the AGN Black Hole Mass Database with the 
NSA and find 19 matches. We require that galaxies 
have photometry in the NSA so we can calculate the 
host galaxy stellar masses consistently with the rest of 
our AGN sample (equation [2). We correct for the AGN 
contribution to the integrated photometry, assuming the 
same power law shape described in section 12.41 How¬ 
ever, for most of the reverberation-mapped AGN, the 
normalization comes directly from L 5100 provided in the 
AGN Black Hole Mass Database. There are a few cases 
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p seudobulges). 

iKormendv fc: Hd l|2013D provide a different way to pre¬ 
dict M/Lk as a function oi B — V color (their equa¬ 
tion 9). This rela tion is based on the m ass-to-light 
ratio calibrations of Unto fc Porting (120131 1 . shifted to 
a dy namically measur e ment zeropoint. Using equation 
9 in IKormendv k, Hd (|2013ll yield s stellar mass e s that 
are systematically higher than the iZibetti et alJ (l2009fl 
masses by 0.33 dex (see Figure [7]). This discrepancy is 
primarily due to different assumed stellar IMFs in the 
m odels and the shift of 0.1 2 58 dex in log M/Ly applied 

( 2009 ) adopt a 
( 2013f l assume 


bv IKormendv fc (1201311. IZibetti et'St 
Chabrier IMF, whereas 


Into fc Portinari 


a Kroupa IMF. The stellar IMF is known to signifi- 
cantly affect the over all normalization of lo g M/Lk (e.g., 
iBell fc de ,Iond [200111 . Here we adopt the IZibetti et al'l 
(j2009f l masses for consistency with the other samples in 
this work. 


Figure 7. Total stellar masses of galaxies with dynamically mea¬ 
sured BH masses derived using iC-band mass-to-light ratios as a 
function oi B — V color from equat ion 9 inIKo rmend v fc Hd Il2ni.'ll l 
versus equation |3] here taken from IZibetti et al.l <20031 . 

where T 5100 is not available, yet we have a measurement 
of broad Ha from the SDSS spectrum. For these, we 
normalize the AGN continuum as described in section 
12.41 Emission lines are not included in the mock AGN 
spectrum for the reverberation-mapped AGN, as these 
measurements are not readily available in many cases. 
We exclude 4 objects (from the initial 19 matches) in 
which the AGN dominates the integrated photometry to 
minimize unreliable stellar mass estimates. BH masses 
and total stellar masses for the 15 reverberation-mapped 
AGNs used in this work are listed in Table |31 

We note that six of the reverberation-mapped AGNs 
in Table [3] are included in our main sample of broad-line 
AGNs with SDSS spectroscopy. As a consistency check, 
we compare the reverberation-mapped BH masses and 
those based on broad Ha emission measured from the 
SDSS spectra (equation 1). For this limited sample of 
six objects, the spectroscopic BH masses are on average 
^ 0.4 dex larger than the reverberation masses. 

3.3. Galaxies with Dynamically Detected BHs 

Benchmark BH masses come from dynamical meth¬ 
ods, which rely on obser vations that spatially re solve the 
BH sphere of influence. IKormendv fc Hoi (j2013f l provide 
an inventory of BH mass measurements based on stellar 
dynamics, ionized gas dynamics, GO molecular gas disk 
dynamics and maser disk dynamics. 

We estimate total stellar masses of galaxies with dy¬ 
namical BH mass measurements using the total abso- 
lute AT-band magnitude s and B — V colors provided by 
IKormendv fc HcJ (1201 3f l in their tables 2 a nd 3 with the 
color- dependent mass-to-light ratio from IZibetti et al.l 
(l2009h given in equation [3] above. The results are 
listed in Table [3 We in clude all objects summarized in 
IKormendv &: Hoi (l2013f) except those with BH mass up¬ 
per limits (2 elliptical galaxies and 2 spiral galaxies with 
pseudobulges), and galaxies without provided B — V col¬ 
ors that are necessary for estimating the galaxy stellar 
masses (2 elliptical galaxiet0 and 3 spiral galaxies with 

Two additional elliptical galaxies do not have B — V colors 


4. LOCAL BH MASS - TOTAL STELLAR MASS 
RELATIONS 

One of the most useful aspects of BH-galaxy scal¬ 
ing relationships is that they provide a way to estimate 
BH mass from more easily measured galaxy properties. 
While the tightest scaling relations appear to be between 
the BH mass and bulge properties in quiescent early-type 
galaxies, we seek to quantify the relationship between BH 
mass and total stellar mass to facilitate work at higher 
redshifts where measuring bulge properties is difficult or 
impossible and BHs are identified via nuclear activity. 

In Figure [5] we plot BH mass versus host galaxy total 
stellar mass for our local sample of 244 broad-line AGN 
and the additional objects described in (jS] A single linear 
relation is disfavored by the data. There is a large range 
in BH mass at a given total stellar mass (e.g., a factor of 
--1000 in Mbh at Mgteiiar ~ 10^°'® Mq). 

Despite a significant amount of scatter in this plot, it is 
clear that at a given total stellar mass, AGN host galax¬ 
ies at z ^ 0 tend to fall below elliptical galaxies and 
spiral/SO galaxies with classical bulges hosting quiescent 
BHs. If the spectroscopic BH masses were shifted down 
(see (I3.2p . this would cause an even larger discrepancy 
between the AGNs and the dynamically detected BHs in 
Figure m Similarly, the AGN host galaxies fall below the 
canonical BH-to-bulge mass relati ons defined by these 
inactive early- type gal a xies (e.g., iHaring fc Rixl 12004 
iMcGonnell fc Ma 120131 : IKormendv &: Hoi l2013f) . BHs 
in galaxies with pseudobulges also preferentially lie be¬ 
low the scaling rela tions based on ellipticals and classi¬ 
cal b ulges (also see lGreene et al.l [20101 : IKormendv et al.l 

[Mil). 

We note that bulge mass and total stellar mass are 
equivalent for elliptical galaxies, which dominate the 
samples used to derive BH-to-bulge mass relations. In¬ 
deed, we find that a linear regression using total stel¬ 
lar mass for the elliptical galaxies and spiral/SO galaxies 
with classical bulges, which have dynamically-measured 
BH masses 1 (I3.3I) . is roughly consistent with standard 

provided bv IKormendv fc Hd Il201.'il ~). however they do have Mtuige 
which is equivalent to total stellar mass since these are ellipticals. 
We include these galaxies in our sample and adopt a total stellar 
mass log Mgtellar = log ll^bulge ~ 0.33, where the offset accounts 
for differences our assumed mass-to-light ratios. 
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Figure 8. Left: Black hole mass versu s tota l hos t galaxy stellar mass. All stellar masses are estimated using color dependent mass-to-light 
ratios presented in IZibetti et al.l (|20Q9I) (see aza and Our sample of 244 broad-line AGN for w hich we estimate vir ial BH masses from 
equation 1 are shown as red points. The 10 broad-line A GN and composite dwarf galax ies from IR eines et al.l (120131 1 are shown as pink 
points (including NGC 4395; Filippenko fc SargenUll989l1 . T he dwarf galaxy RGG 118 iR eines et^L n^olSI l hosting a ~50,000 M (^ BH 
iBaldassare et al1l2015ll is the dark green point, and Pox 52 IBarthetjulgOT^Tl^ornton et‘‘ani200^ is the light gr een point (see jXTJ. 
Fifteen reverberation-mapped AGN with BH masses taken f rom IBentz Ik Kat3 II2015I 1 are shown as purple points (see i|3.2ll . Dynamical BH 
mass measurements are taken from IKormendv &: Hd (12013 1 and shown as blue (elliptical galaxies), turquoise (S/SO galaxies with classical 
bulges) and orange (S/SO galaxies with pseudobulges) points. The gray error bar indicates uncertainties in stellar masses for all points, and 
single-epoch spectroscopic BH masses. The gray lines show various Mbh vs. M^uige relations based on ellipticals and spiral bulges with 
dynamical BH mass measurements. The lKormendv Hi Hoi (12013 1 “scaled” relation has bulge masses scaled down by 0.33 dex to account 
for differences in our assumed mass-to-light ratios (see ij3.3ll . 


bulge mass relations, albeit with more scatter (see be¬ 
low). 

Thus, it appears that a separation exists between our 
sample of uniformly selected AGN hosts (Sj2]), and el¬ 
lipticals and classical bulges. We anticipate that using 
or extrapolating the canonical BH-to-bulge mass scaling 
relations to interpret samples of galaxies with uncertain 
morphological classification, or AGN hosts, may lead to 
erroneous inferences. 


4.1. The BH-to-Total Stellar Mass Relation for Local 

AGNs 

We plot log Mbh versus log Mgteiiar for the AGNs 
alone in the left panel of Figure El We first use a non- 
parametric method to help visualize the data and demon¬ 
strate that there is indeed a correlation between BH mass 
and total stellar mass for local A GNs. We use the kernel 
density estimation technique (e.g.. lSilvermanl[T98^ to es¬ 
timate the density function in the log Mgteiiar — log Mbh 


plane from the observed data for all AGN13- Each data 
point is represented by a two-dimensional normalized 
Gaussian kernel. The smoothing parameter (e.g., a for 
a Gaussian) is set to 0.3 and 0.5 for log Mgteiiar and 
log Mbh, respectively, and reflects the measurement un¬ 
certainties for the majority of our sample (where masses 
are in units of Mq). The individual kernels are then 
summed to produce the kernel density estimate (left 
panel of Figure ED- The kernel density estimate is subse¬ 
quently normalized for each log Mgteiiar independently to 
construct the conditional probability distribution func¬ 
tion (PDF), p(logMBH|logMsteUar), which illustrates the 
dependence of BH mass on total stellar mass for our sam¬ 
ple of AGNs. The right panel of Figure El shows the re¬ 
sulting PDF, where the lines correspond to the median 
and standard deviation as a function of log Msteiiar- 

This non-parametric method nicely illustrates a cor- 

® For individual AGN with multiple BH mass estimates, we 
include only one data point with priority given to reverberation 
masses when available (e.g., NGC 4395). 
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Figure 9. Left: Black hole mass versus total host galaxy stellar mass for local AGNs with the kernel density estimate (see text) shown in 
grayscale. Contour levels are at (1/2)’^ x the peak value, where n = 1 to 5. Right: Conditional PDF p(logMBH l^og^stellar) computed by 
normalizing the kernel density estimate at each log Mgteiiar- The middle line indicates the median of the PDF as a function of log Mgteiiar 
and the outer white lines show the standard deviation. 


relation between log Mbh and log Mgteiiar for our sam¬ 
ple of AGNs. However, the shape/slope of the relation 
may be different for the population of local AGNs since 
the conditional PDF is based on data that is suscepti¬ 
ble to selection biases that are particularly severe at low 
masses (see gSH). While we should not immediately as¬ 
sume these data are well described by a linear relation, 
the sample Pearson correlation coefficient indicates that 
a linear relationship between log Mbh and log Mgteiiar 
is a reasonable description of the data; r = 0.54 with a 
probability p < 10 “® that no linear correlation is present. 
We therefore use a line to parameterize the AGN rela¬ 
tion. 

We take a B ayesian appro ach to linear regression using 
the method of iKellvl (l2007rFI . which accounts for uncer¬ 
tainties in both log Mbh and log Mgteiiar- To facilitate 
comparison with other studies, we parameterize the re¬ 
lation as 

l 0 g(MBH/MQ) =a + 13 log(Mgtellar/lO^^M 0 ) (4) 
and find 

a = 7.45 ±0.08;/?= 1.05 ±0.11. (5) 

The quoted slope and intercept are given by the median 
of 10,000 draws from the posterior probability distribu¬ 
tion of the parameters. The errors on the linear coef¬ 
ficients are correlated and the reported values are de¬ 
termined from the Icr error ellipse. The rms deviation 
of the BH mass measurements from the relation is 0.55 
dex, and incorporates both our adopted measurement 
errors of 0.50 dex and a best-fit intrinsic scatter of 0.24 
dex (added in quadrature). The intrinsic scatter may 
be larger if our measurement errors are overestimated. 
The linear relation for the AGN host galaxies is shown 
in the left panel of Figure [101 We note that our results 

® linmix.err.pro in the IDL Astronomy User’s Library. 


do not change significantly when using only our primary 
sample of broad-line AGNs ([2]) a nd e xcl udin g the addi¬ 
tional AGNs described in sections l3Tl and IT^ The slope 
and intercept derived from our uniformly selected sample 
agree with those in equation within the Itr uncertain¬ 
ties. 

The BH-to-total stellar mass relation for the ellipticals 
and classical bulges is shown in the right panel of Fig¬ 
ure [10] for comparison. The corresponding coefficients in 
equation [4] are given by 

a = 8.95 ±0.09;^ = 1.40 ±0.21, ( 6 ) 

and the intrinsic scatter is 0.47 dex (rms = 0.60 dex). 

4.2. Implications for BH-to-Stellar Mass Fractions 

Figure [TT] shows BH mass fractions as a function of 
stellar mass for our two total stellar mass relations at 
z ^ 0 (AGNs and dynamically measured BHs; equa¬ 
tions HID, as well as some standard BH-to-&u/qe mass 
relations le.g.. iKo rmendy fc Holl20T^ iMcGonnell fc Mai 
1201,81 : iHaring fc Rixll2004l l. The BH-to-total stellar mass 
fraction given by the AGN relation is Mbh/M gteiiar 
0.02% to 0.03% across the stellar mass range 10® < 
Msteiiar/M 0 < 10^^. This is markedly smaller, by 
roughly an order of magnitude, than the BH-to-bulge 
mass fractions derived from quiescent early-type galax¬ 
ies that are commonly used as references. 

For instance, the BH mass fraction given by our AGN 
total stellar mass relation is ^ 19 to 56 times lower than 
th at given by the canonic al BH-to-bulge mass relation 
of iKormendv &: H'3 (|2013r ) across the stellar mass range 
10 ® < Msteiiar/M 0 < 10 ^^, accounting for differences in 
our assumed stellar mass-to-light ratios (see >13.31) . Using 
total stellar mass rather than bulge mass for the host 
galaxies of dynamically measured BHs in ellipticals and 
classical bulges also results in a BH mass fraction that 
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log (M3teii3/Mg) log (M3,g||3/Mg) 

Figure 10. Left: The BH-to-total stellar mass relation for local AGNs (dark red line; equations [4] and [Sj . The light gray shaded region 
accounts for the errors in the slope and intercept of the relation and the dark red dotted lines indicate the rms scatter of points around the 
relation (0.55 dex). The gray error bar indicates uncertainties in stellar masses for all points, and single-epoch spectroscopic BH masses. 
BH mass errors for the reverberation-mapped AGN are shown on the (purple) data points. Right: Same as the left panel, but for the 
inactive sample of elliptical galaxies and spri al/SO galaxies with classical bulges (equations |4] and |6} . The dark blue line indicates our 
relation derived using total stellar mass f 113.311 . 


is roughly an order of magnitude larger than that of the 
AGN host galaxies. We thus urge extreme caution when 
using the canonical BH-to-bulge mass scaling relations 
as a proxy for BH-to-total stellar masses since this may 
lead to a biased interpretation. 

4.3. Systematic Uncertainties 

Based on the z ~ 0 samples presented in this work, 
we have shown that AGNs occupy a region in the 
log Msteiiar — log Mbh plane below that populated by 
ellipticals and classical bulges. Here we consider sys¬ 
tematic uncertainties in stellar and BH masses that may 
affect this empirical result. We note that differences in 
distance scales are negligible. The range of assumed Hq 
across all samples represented in Figure [5] varies from 70 
to 71 km s~^ Mpc“^. 

We do not expect that systematic uncertainties in stel¬ 
lar masses between the samples shown in Figure [10] can 
account for the two different relations. First, all stellar 
masses have been estimated in the most consistent man¬ 
ner feasible. That is, we use color depend ent mass-to- 
light ratios provided by iZibetti et al.l (|2009D with either 
a combination of g and i-band SDSS data (all AGN ex¬ 
cept Pox 52), or B, V, and AT-band data (galaxies with 
dynamical BH masses and Pox 52). As described in H2.4\ 
we have accounted for any AGN contribution when calcu¬ 
lating stellar masses so we do not think AGN host galaxy 
masses are significantly overestimated. Moreover, in or¬ 
der to bring the sample of AGNs onto the upper relation 
by shifting their stellar masses, the stellar masses would 
need to be reduced by more than an order of magnitude 
(see Figure [5]). 

The virial BH masses estimated for our sample of 
broad-line AGNs are quite indirect and subject to vari¬ 
ous uncertainties. For example, the broad-line region ge- 



Figure 11. BH mass fractions (given as a percentage of the stellar 
mass) as a function of stellar mass. Our local AGN relation, where 
^stellar = ^total is shown as a dark red line. Our total stellar 
mass relation for quiescent BHs in early-type galaxies is shown as 
a dark blue line. Bulge mass relations from the literature are shown 
as gray lines. 

ometry and orientation certainly varies between objects 

(lKollatschn^l200^lBentz et al.ll2009l : iDennev et al.ll201^ 

iBarth et al.ll2011l) . yet we apply a single geometric scal¬ 
ing factor since we do not have this information for the 
individual broad-line AGNs in our sample. There is also 
the possibility that there are nongravitati onal contribu - 
tions to the measured gas velocities (e.g.. lKrolikl[20nl . 
although this would lead to systematically overestimated 
BH masses. 




















12 


Reines et al. 


The lower relation defined by the broad-line AGNs 
has a normalization that is ~1.2 dex lower than the 
upper relation at Msteiiar = 10 ^°Mq. Across our 
sample, uncertainties in the virial BH mass estimates 
are expected to be on the order of ^ 0.5 dex (e.g., 
IVestergaard fc Petersonll^OObl : LShenll201,‘lfl . which is con¬ 
siderably less than the offset in BH mass between the 
two relations. Any reasonable variation in the virial fac¬ 
tor can also be ruled out as producing artificially low BH 
masses for the broad-line AGNs. In order to get AGN BH 
masses to fall on the upper relation by changing the virial 
factor, < / > would need to be > 40 (e > 10 in equation 
1). Finally, we note that dynamically measured BHs in 
galaxies with pseudobulges, as well as the reverberation- 
mapped AGN, overlap our sample of broad-line AGN. 
For all of these reasons, we conclude that uncertainties 
in virial BH masses alone are not artificially producing a 
lower relation in the log Mgteiiar — log Mbh plane for the 
AGN. 

We note that we have discarded 9 objects from our 
sample of AGNs in which the luminosity of the AGN 
dominates the total integrated photometry (host -|- 
AGN). Our motivation for this was to minimize unre¬ 
liable stellar mass estimates. For a fixed stellar mass 
and Eddington ratio, this could bias us against large 
Afen/Alsteiiar- However, given that only ~ 3% of the 
AGNs were removed from our sample, we do not think 
this has impacted our results in any appreciable way. 

4.4. Possible Origins for the Separation Between AGN 
Hosts and Ellipticals/Classical Bulges 

We now turn our attention to possible origins for 
two BH-to-total stellar mass relations: one comprising 
AGNs, the other ellipticals and classical bulges, with 
pseudobulges predominantly overlapping the AGNs. It is 
worth noting that at least some of the galaxies with pseu¬ 
dobulges are active, as they are selected as being masers. 
Others, like the Milky Way, are obviously inactive. 

Given that the lower relation is defined by AGNs and 
the upper relation is defined by quiescent BHs with dy¬ 
namical BH mass measurements, it is reasonable to con¬ 
sider if nuclear activity (or lack thereof) may be par¬ 
tially responsible for the existence of two separate re- 
altions. On the one hand, there are reasons to ques¬ 
tion the importance of nuclear activity since the samples 
defining the two relations are each fraught with their 
own selection biases. Dynamical BH mass measurements 
are severely biased towards nearby, massive and dense 
galaxies wh ere the BH sphere of influ e nce can be re¬ 
solved (e.g.. Ivan den Bosch et al.ll2015ll . iGiiltekin et al.l 
(IMHI showed that dynamical BH masses are not bi¬ 
ased high in the Mbh — o'* plane for very large galaxies 
(< CT* 268 km s“^), however we are probing signif¬ 
icantly lower masses. There is no reason that galaxies 
with quiescent BHs should not overlap with the AGN 
host galaxies in Figure[5] we just cannot detect such BHs. 
One the other hand, we see in Figure[8]that we do not de¬ 
tect AGNs with BH masses as large as the quiescent BHs 
at a given stellar mass, suggesting nuclear activity is an 
important factor. If there were AGNs with such large BH 
masses, we should see them since they would be bright. 
Apparently, the larger BHs (at a g iven stellar mass) are 
not shining as AGNs (see also, e.g.. lHeckman et ^12004 
lMerlonill2004[l . 



log (Msteiia/Mg) 

Figure 12. Log Mbh versus log Mgteiiar with p oints color-coded 
by approximate morphological type (see text in i|4.4l l. AGNs are 
plotted as stars and dynamically detected BHs are plotted as cir¬ 
cles. The dark red line indicates our BH-to-total stellar mass re¬ 
lation based on local AGNs (equations and O. The dark blue 
line indicates our BH-to-total stellar mass relation based on ellipti¬ 
cals and galaxies with classical bulges hosting dynamically detected 
BHs (equations 13 and |6]l . 

Another potential factor is differences in host galaxy 
properties. A comparison of the Hubble types of 
the galaxies with dy namically detected BHs (from 
iKormendv fc Holl2013l) supports this notion. For exam¬ 
ple, the classical bulges are mostly hosted by early-type 
SO galaxies (16/20 have Hubble types of SO or SBO), and 
they overlap with the elliptical galaxies. Alternatively, 
the pseudobulges that tend to overlap with the AGN 
host galaxies are more commonly found in later-type spi¬ 
ral galaxies (2/17 are found in SBO galaxies, the rest are 
in spirals). 

We find that a significant fraction of local broad¬ 
line AGN host galaxies are also spirals/disks, simi¬ 
lar to what has previously been found in studies of 
moderate-l uminosity AGN ho s t galaxies out to z ^ 
3 (e.g., iGabor et al.l 120091: iScha-winski et al.l 120111 : 
iKocevski et al.l 120121 : iBennert et al.l I2015D ~ We ob¬ 
tained approximate Hubble types for a subset of the 
AGN host galaxi es using automated morphologi cal clas¬ 
sifications from [Huertas-Gompanv et al.l (|2011f) . The 
automated classifications from this work have been 
shown to correlate well with visual classifications. 
Rath er than assigning a singl e morp hological classifica¬ 
tion, iHuertas-Gompanv et al.l (1201 Ih provide probabili¬ 
ties of being in four morphological classes (E, SO, Sab, 
Scd). For the AGN host gala xies with matches in the 
[Huertas-Gompanv et al.l (1201 Ih catalog (129 objects), we 
simply adopt the class with the highest probability. 

In Figure [12] we plot log Mbh versus log Mgteiiar with 
points color-coded by morphological class. We also in¬ 
clude the galaxies with dyna mical BH masses , convert¬ 
ing Hubble types provided bv IKormendv fc Hoi (|2013D to 
either E, SO, Sab, or Scd. Based on these rough morpho¬ 
logical classifications, it is clear that a significantly higher 
fraction of spiral galaxies lies close to the lower relation 













































Local Mbh — Af* Relations 


13 


compared to the upper relation. These spiral galaxies 
will have less prominent bulges than the galaxies on the 
upper relation, which are mostly bulge-dominated ellip¬ 
tical and SO galaxies. 

The AGN host galaxies on the lower BH-to-total stel¬ 
lar mass relation could conceivably follow the canonical 
BH-to-bulge mass relatiorQ if the classical bulge masses 
for the AGN hosts were on average only ^ 5% of the 
total stellar masses. T his echoes the re sult of the model 
for bulge evolution bv iLu fc Mol (j2015r i. where their ref- 
erence model ( Mode l III) is a very good fit to our data. 
iSanghvi et al.l (j2014l l also find that low-mass quasars lie 
below the extrapolation of the local BH-to-bulge mass 
relation, but with a correction for the disk they obey it. 
iGaplar et al.l (|2014ll advocate a lower MBn/A^steiiar ra¬ 
tio for local active galaxies compared to inactive galaxies 
that have “quenched” at earlier times. 

5. CONCLUSIONS AND DISCUSSION 

In this paper we have studied the relation between 
BH mass and total stellar mass for nearby galaxies 
(z < 0.055), including both galaxies with inactive BHs, 
with dynamical BH mass measurements, and galaxies 
with an active BH, with mass measurements based on 
reverberation mapping or single-epoch virial estimates. 
Inclusion of the latter sources allows us to extend our 
sample to low BH masses, and also to use the same tech¬ 
nique used at higher redshift. Moreover, our stellar mass 
measurements rely on mass-to-light ratios, as is routinely 
the case for higher redshift samples. Therefore we build 
a local analog of higher redshift samples, but where we 
have better c ontrol on systematic s. Our work is com¬ 
plementary to iLasker et al.l (j2014fl , where they measure 
in a uniform way the total luminosities of a sample of 
galaxies with dynamically-measured BHs, in that we try 
to provide a benchmark for observational or theoretical 
studies where detailed inform a tion o n bulge properties is 
not available. iBennert et al.l (j2015f l also took a similar 
approach towards the correlation with velocity disper¬ 
sion. One important caveat in our analysis is that the 
AGNs we have studied are moderate-luminosity Seyferts 
(41.5 < log Lboi ^ 44.4). At high-redshift, statistical 
samples are often biased towards more luminous AGNs 
and quasars. 

In Figure m we plot BH mass versus total stellar mass 
for galaxies with 10®'® < Afgteiiar ^ 10^^ Mq. A single 
linear fit to the data is disfavored. Rather, the AGN 
host galaxies define a relation that has a similar slope 
(Mbh oc Mgteiiar) to early-type galaxies with quiescent 
BHs, but a normalization that is more than an order 
of magnitude lower (Figures [10] and [H]). The differ¬ 
ent normalizations may be partially due to active ver¬ 
sus inactive BHs, but can also be attributed to differ¬ 
ences in host galaxy (bulge) properties as discussed in 
1 14.41 We caution that using the z = 0 benchmark BH-to- 
bulge mass relations for AGN host galaxies, or assuming 
Mbuige = Mtotai, may lead to severely biased interpre- 
tations. For example, [Graham fc ScottI (1201511 assume 
bulge mass equals total stellar mass for t he 10 broad- 
line AGN and composite dwarf galaxies in iReines et al.l 

^ We have scale d the BH-to-bulge mass relation of 
IKormendv fc HH II201,'1I1 to account f or ou r different assumed stellar 
mass-to-light ratios as discussed in ED 


(|2013ll . which in part leads them to conclude that AGN 
host galaxies follow a steeper, quadratic-like relation be¬ 
tween BH mass and bulge mass. 

Our work also has important implications for cosmo¬ 
logical simulations that are tied to the local BH-to-bulge 
mass relations. In most cosmological simulations that 
produce statistical samples of BHs and AGNs, i.e., large 
un iform-volume simulation s (e.g., MassiveBlack I and 
H, iDi Matteo et al.l (1201211 : iDeGraf et al.l (1201411; IHus - 
tris. iSiiacki et al.l (12015); Eagle, ISchave et al.l (j2015D : 
Horizon-AGN. jPubois et fo] (j2014[ B the resolution is lim¬ 
ited to ^kpc scales, making bulge-disk decomposition 
for low-mass galaxies unreliable. The typical approach, 
in fact, is to estimate th e total stellar mas s within 
the st ellar half-mass radiu s (ISiiacki et al.ll201^ or twice 
that (jPeGraf et al.l l2014|i , or extrapolate a fit to the 
mass profile of the bulge inferred from kinematic data 
(jSchave et al.ll201^ . If the simulations do not select el¬ 
lipticals or galaxies with classical bulges and perform a 
bulge-disk decomposition, using the relat ions published 
in th e literature for BH - bulge mass (e.g.. iHarin g fc Ria 
2f)nl IMarconi fc Huntl I^M IMcConnell fc Mai 1201,31: 


Kormendv fc Hcill2013f i as a benchmark for the compari¬ 


son between simulations and observations at z = 0 (and 
perhaps beyond) would not be appropriate. Our method 
is easier to implement in the analysis of simulations, and 
can help disentangle issues related to BH growth and 
AGN properties. 

Whether or not the BH-to-total stellar mass relation 
extends to even sma ller masses bears directl y to th e ori- 
gin of BH ‘seeds’. iVolonteri fc NataraianI (1200911 and 
Ivan Wassenhove et al.l (1201011 suggested that if BH seeds 
are massive, e.g., 10® Mq, as predicted by ‘direct col¬ 
lapse’ models, the low-mass end of the relation between 
BHs and galaxies (they specifically referred to the corre¬ 
lation with the velocity dispersion, but the result would 
hold for the stellar mass as well) flattens towards an 
asymptotic value, creating a characteristic ‘plume’ of un¬ 
grown BHs. Vice versa, if BH seeds are small, e.g., 10^ 
Mq, as predicted by models related to the first genera¬ 
tion of stars, the expectation is that the observable scal¬ 
ing laws would not see the asymptotic value (the ‘plume’) 
because it lies at masses below those that can be probed 
observationally. 

Finally, extending the sample at stellar masses ^ 
10® Mq is fundamental to interpret results for much 
higher redshift galaxies with similar masses. Searches 
for AGN in galaxies with stellar masses ^ 10® Mg) at 


z > 6 have found very few, if any, BHs (iWillott 


Fior£jt^L|j2012yCowj£e1^jljj20f2 et al 


2011 


2013 


Giallongo et al.ll2015l: IWeigel et al.ll2015|i . If our sample 

is a good representation of the local universe, and the 
same relation based on local AGNs holds at high-redshift, 
it explains why we do not easily detect BHs at high-z: 
their masses would lie below the extrapolation of the lo¬ 
cal BH-bulg e mass relation wh i ch is normally used as a 
benchmark. IVolonteri fc Stard (|2011h suggested that ob¬ 
servations could indeed be explained with a BH-stellar 
mass correlation eithe r steeper than at z = 0, or with 
a lower normalization. iDubois et al.l (1201511 suggest that 
the growth of BHs in low-mass galaxies (galaxy mass 
< 10^° Mq and bulge mass < 10® Mq) is stunted be¬ 
cause of supernova feedback, which hinders accumulation 
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of gas in the nucleus until the potential well of the bulge 
and galaxy become deep enough. In forthcoming work, 
we will explore the consequences of our results on the 
interpretation of high-z BH populations, and the link to 
their hosts. 
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